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Abstract 

We exploit the variation of the atomic interaction in order to move ultra-cold atoms across 
an AC-driven periodic lattice. By breaking relevant symmetries, a gathering of atoms is achieved. 
Accurate control of the gathered atoms' positions can be demonstrated via the control of the atomic 
localization process. The localization process is analyzed with the help of the nonlinear Floquet 
states where the Landau-Zener tunneling between states is observed and controlled. Transport 
effects in the presence of disorder are discussed. 

PACS numbers: 03.75.Lm, 05.60.-k, 63.20.Pw 
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I. INTRODUCTION 



In nonlinear discrete lattices, spatial discreteness and nonlinearity constitute the two 
main ingredients necessary for the existence of localized excitations. The existence and 
dynamics of these excitations can be understood from the analysis of their energy spectrum. 
In this spectrum, the localized states lie in gaps between bands of extended states. 

The spectrum, and therefore the localized states, depend very much on the strength of the 
nonlinearity since the position and width of the bands and gaps vary with this parameter. 
In this respect, the nonlinearity strength can be used as a control parameter for tuning 
localized excitations. An ideal experimental testing ground for the application of this type 
of control is that of a Bose-Einstein condensate (BEC) trapped in a deep optical lattice, 
where much of the nonlinear and discrete effects can be probed in optical lattices jlj]. 

Further control of matter waves (BEC) in optical lattices has been demonstrated with 
the use of fields. AC-fields can modify the properties of the matter waves, by adjusting 
the parameters of the fields This is an intense and growing area of research, where 
examples of striking phenomena, such as transitions of a superfluid to a Mott-Insulator [2I 
and generation of directed transport (4, 5], have found an experimental footing jf], Q]. 

In this respect, a combination of AC fields together with management of nonlinearity 
opens new avenues for exploration of the new phenomena in BECs trapped in optical lattices. 

In this work, we consider a general management of the atomic interaction for a BEC in 
presence of an AC driving field. Manipulation of the atomic interaction is usually referred 
to "management of Feschbach resonances" |{|. At these resonances, the scattering length 
undergoes large fluctuations with values going from positive to negative as the magnetic 
fields are tuned J9], thus opening many possibilities for the control of the atomic interaction. 
"Management" here refers to manipulation of DC and AC parts of the atomic interaction. 
The AC part is locally applied to induce motion of the atoms towards the perturbed sites. In 
the light of new experiments, local variation of the atomic interaction is feasible by spatially 
modulating the interaction strength on a short-length scale [lo|]. Here, rather than focusing 
on the localization process itself, we analyze how to get transport of particles out of the 
localization process; i.e., how to transport cold atoms within optical lattices by controlling 
the location of the localization sites. 
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II. MODEL 



Consider a lattice with periodic boundary conditions under the action of an AC force. 
The dynamics can be described by the equation 

i^- = E n ip n + C(i/j n+1 + V>n-i) +gip n N n - il) n F n {t), (1) 

where E n is the energy at the site n. N n = N\ip n \ 2 is the population of atoms at the i th 
site. N is the total number of atoms, that for convenience is normalized to 1. For the sake 
of dimensionless units, we set h = 1. g accounts for the interaction strength between atoms, 
which is a function of the s-wave scattering length. This in turn can be tuned by means of 
magnetic fields [9|. More recently, control of the scattering length has been realized with 



a laser |lOj, allowing more accuracy and a less significant loss of atoms. F n (t) in EqJT] is 
a periodic function that mimics a flashing potential |6| for a discrete lattice. This can be 
written as 

F n (t) = (-l) n f(t) = (-l) n hajx{tJt) i (2) 

where u is the driving frequency and h is the amplitude. The term (—1)™, accounts for the 
spatial periodicity of the potential. Eq.((T]) with fl2]) satisfies periodic boundary conditions. 

In Eq.([T]) the atomic interaction strength g can include CD and AC terms whose manip- 
ulation is usually referred to as "management of Feschbach resonances (MFR)" |sj]. MFR 
las proven to be an effective tool, on the mean field level, for manipulation of cold atoms 



. It has a 



of atoms rjj . 



so been suggested as way to control the tunneling process of a discrete number 



Notice, moreover, that the control of atomic population imbalances in a double well poten- 
tial has been realized by varying parameters in the system including the atomic interaction 



strength between atoms 12j. A similar kind of control has been proposed for controlling 
the wavepacket spreading in a disordered lattice [kJ • Upon varying the ramping speed of a 
linear time-dependent strength of the nonlinear term, the authors find optimum localization 
in the lattice in presence of disorder. 

From the perspective of a cold-atom system, localization, (namely gathering of atoms at 
few sites) implies a motion of atoms across the lattice, meaning a transport of particles. 
Notice that the process here is different from that of moving solitons, where the motion 



is achieved by kicking localized matter waves 



14|. Thus, controlling the location of the 
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FIG. 1. (Color online) (a) Quasienergy e vs. go for the two modes system [Eqs. (jH) and ©]. 
Green and red solid curves: U = 0; black dashed curve: JJ max = 2. Upper inset: Examples of 
pulse shape: U vs. go- two-slope (solid line) and Sine-square (dashed line). Lower inset: Average 
of population imbalance (S) vs. go for JJ max = 2. (b) (S) vs. phase 9, JJ max = 2. The average 
(...) is realized over a period T = 2ir/oj. Computations are performed with the sine-square shape 
depicted in the upper inset of (a). The parameters are C = 1, u = 2ir, h = 1, E = 0, go = 3. 

localization sites allows us to control the transport of the atoms. With this aim, we consider 
DC and AC terms for g, namely 



where 9 is a phase. It is expected that the periodic field in Eq.(j2]), along with management 
of the atomic interaction, may help to control localization of atoms at specific sites of the 
lattice. 

A. Two-mode system 

To gain an insight into the dynamics, it is convenient first to analyze the dimer lattice, 
which can be written as 




(3) 
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Et/> 1 + Cik + 94>iNi + il>if{t), 
-E^ 2 + C^ x + g^N 2 -f 2 f{t). 



(4) 
(5) 



Here the energies are set as E% = E and Ei = —E, where E can be seen as a bias. 
Eqs. (jl]) and (j5]) with E = f (t ) = have been used to describe the self-trapping transition 



a critical value g c 



of two weakly cou pled BEC 15], This phenomenon appears when the nonlinearity exceeds 
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16|. 



To understand this phenomenon in the presence of an AC field fit), it is convenient to 
study the eigenfunctions of the extended Hilbert space of T-periodic functions, similar to 
that implemented in the linear regime (see Ref. 13]). Conventional Floquet theory fails in 
the presence of a nonlinear dependence on the wave functions. Nonetheless, Floquet states 
of the linear system can be traced into the nonlinear domain and new periodic orbits are 
found, which are then called Nonlinear Floquet states 12j, ll8H20| . The periodic solutions in 
the linear regime g = fulfill the relation \ipp(t)} = e~ i£ ^ \<f>p(t)), where \<f> p (t+T)) = \<f>p(t)) 
and ep is the quasienergy value. 

To characterize the localization of atoms in a two-mode system, we use the quantity 
5 = \ipi\ 2 — \tp2 \ 2 i which is the atom population imbalance between the two sites. 

Equations (jl]) and (jSJ) have two Floquet states for g = with zero population imbalances, 
namely 5 = 0. When increasing g above some g c , one of the two states bifurcates into 
three new states, as shown in FigOJi. Two of these new states are degenerate, with opposite 
population imbalances. These are impossible to reach by continuation from the linear regime, 
because of the breaking of the adiabatic condition at the branching point 3]. Let us now 
analyze the symmetries of Eqs. (jl]) and (jSJ) in the presence of MFR (see Eq.(j3])). 

If U = in Eq.(j3]), Eqs. (jl]) and (jSJ) become invariant under the transformations: 



and 



5x : Vi t^t + T/2, 



S2 '■ 4>i — > t — > —t, complex conjugation, 



(6) 



(7) 



where 5i corresponds to a generalized permutation symmetry. 

On the contrary, for U 7^ 0, S± is always broken, and 52 is broken if 9 7^ ±7r/2. 



5 



0.5 

<S> 
-0.5 

-0.2 -0.1 0.1 0.2 
E 

FIG. 2. (Color online) Average population imbalance (S) vs. bias E for U max = 2, h = 1 (black), 
h = 2 (red); go = 4, 9 = 0. We compute in the limit a = 0. Computations are performed with 
the sine-square shape depicted in the upper inset of FigQJa). Inset: Quasienergy e vs. go- The 
parameters are C = 1, h = 1, E = 0.1, U = 0. 

Thus, taking Eq.(j3]) in Eqs. (jlj) and ([S]) with 9 ^ tt/2, breaks both Si and £2 symmetries. 

To follow states from the linear regime with zero population imbalances, it is convenient 
to modulate the AC nonlinear term U in Eq. ([3]) as a function of go, namely U(go)- In the 
following, we assume U(go) has a pulse shape. As a result of this modulation, a new energy 
branch is created off the branching point that continuously join states with zero or little 
population imbalance to those states of strong population imbalances in Fig{TK The pulse 
shape appears not to be very important, since the plots of the quasienergy curves obtained 
for pulses with a two slope- and sine-square shapes (see upper inset of FigfTJi) show no 
significant difference. From here onwards, we take on the two slope-shape. To reach those 
states with strong population imbalances, we simply ramp go linearly in time, viz. 

9o{t) = at, (8) 

where a is the speed of variation or ramping rate. 

Interestingly, changes of the phase 9 creates states with opposite population imbalances, 
as shown in lower inset of Figflji. A full scan of the phase shows that the population 
imbalance becomes zero for 9 = ±n/2 only (see Fig. [lb). This is in agreement with the 
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symmetry analysis from Eqs. flS]) and (J7J). States with opposite population imbalances can 
therefore be selectively targeted when slowly changing go. Moreover, the smoothness of the 
new energy branch in the vicinity to the bifurcation point allows a significant increase in 
the ramping speed a. This is because a sharp variation in quasienergy is usually interpreted 
as a strong interaction with a very near eigenstate 23j. If the energy variation is smooth, 
high a values can be taken into account. However, for too high a, the system may hop from 
one state to other states with zero or little population imbalance. This implies a loss of 
localization. This jump brings to mind the Landau-Zener tunneling observed in the linear 
regime 23] . 

Furthermore, Eqs. (j3J) and §J5§ can describe the dynamics of opposite momenta for a con- 
densate [20J]. A generalization to an extended lattice, together with Landau-Zener tunneling 
between states of different momenta, may significantly change the transport of the atoms 

The above technique comprises a rather general MFR (EqlH]) for the atomic interaction. 
Thus, we can partly summarize the main effects on the spectrum of quasienergies when 
the above generalized management of the atomic interaction is applied: first, it creates a 
new quasienergy path that circumvents the branching point, second, the new path allows a 
rapid passage, and third, positive or negative strong population imbalances is reached upon 
choosing the phase 9. 

So far, the analysis has been focused for a zero bias, i.e. E = in Eqs.(H| and (jHJ). 
For a nonzero bias E 7^ 0, degeneracies are lifted as shown in inset of FigJ21 In this case, 
localization can also be achieved by slowly changing go in time. In the same spirit as the 
functioning of a ratchet system 24|, here we are interested in how localization is induced 
by means of AC fields and, more importantly, how this works against a bias. Our results 
show that the induced localization appears to work against a bias as shown in Figj2j An 
averaged nonzero population imbalance is induced for E = 0, whose effect survives for finite 
values of E. The offset is shifted to the right as the amplitude of the field f(t), h, increases. 
Likewise, one can induce the same effect to the left by changing 9 in f(t), Eqs.(@| and (E]). 
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FIG. 3. (Color online) Nonlinearity strength go vs. quasienergy e for n = 4. h = 4. Dashed line is 
a quasienergy state that results from continuing the linear Floquet state into the nonlinear domain, 
in the presence of an AC pulse of the nonlinearity strength, with the profile shown in the inset. 
Inset: Pulse profile U vs. go, uj = 4ir, C = 1. 



B. Lattice 



As mentioned above, states bifurcate in new branches, representing new solutions, as 
the nonlinearity strength increases. Hence, in order to describe the many possible solutions 
that appear in the system of Eqs. ([T|) and (J2J), as the nonlinearity increases, a shorthand 
description for each solution of the branch is convenient. Such analysis is usually realized 
in the limit of q — > oo, where one can estimate the asymptotic solution -0 = (ipi, t/> 2 , ...) for 
every branch 25 1. 

Rather than determine the solution ift = (^i, ip2, ■■■), here we are interested in the 
modulus-squared of the wave function in the sites, i.e, \ipi\ 2 . That is, the number of atoms 
localized at each site i. Thus, we take a modified shorthand representation, where now "•" 
accounts for the occupied sites with a 1/K fraction and "o" for a zero fraction of atoms of 
the corresponding state in the limit (7—7-00. Here K is the number of nonzero ipi. 

In the following, we consider a four- mode lattice. In this case, EqfTJfor four sites remains 
invariant under permutations of the indices 1 — > 3 and 2 — > 4. Likewise, symmetries given by 
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FIG. 4. (Color online) Lower: Sketch of the ramping process. Upper: Distribution of atoms in a 
ring of connected sites at the beginning and end of every ramping process. (For meaning of colors 
read text). 

Eqs. flSD and ([7]) are also applicable. Figure [3] shows a section of the quasienergy spectrum vs. 
nonlinearity strength. Of particular interest are the states {• • • •} with equal distribution 
of atoms across the lattice and {• o o o}, with all the atoms localized in one site. 

We consider now the MFR [Eq.([5])] which creates as mentioned above, a new branch that 
join the zero-population-state with that of very high population imbalance as go increases 
(Fig 131 dashed line). 

Interestingly, as the pulse of the AC nonlinearity strength U(go) is applied on two sites 
only, the number of atoms tend to localize in one of these two sites when ramping go- 
Meaning that a motion of atoms from all the other sites to this new populated site has 
taken place. This can be fully controlled with the phase 9 [see Eq.([3])]. In this respect, 
the combined action of ramping go along with an AC perturbation of g [see Eq.([3])] can, in 
principle, help to control the movement of atoms along the lattice. Figure 0] schematically 
shows a method to move atoms from one lattice site to another. The lower graphic depicts 
two stages in time, each one consisting of two ramping process of go, which are described as 
follows: As time proceeds go increases until it reaches a maximum value g™ ax at r = g nax /a, 
where a localization of atoms takes place. Afterwards go is ramped down in time until go = 
at 2r (at this point the number of atoms are equally distributed). At the second stage, the 
AC perturbation is applied on two sites, different from the prior ones, thus making atoms 
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gather in a different new lattice location. In doing so, we have effectively achieved moving 
atoms from one lattice site to another at time 3r. The entire process of moving atoms 
to different sites in a circular ring is depicted in the upper part of FigJH by a sequence of 
drawings. Each drawing represents a configuration of localization of atoms in the ring at the 
end of every ramping process. The red color on one site stands for a large concentration of 
atoms at this site, whereas absence of color represents a zero or very low number of atoms. 
Drawings with brown color at every site means an equal distribution of atoms. Notice that 
this process can be repeated over and over, to shift atoms to any target site of the lattice. 

A remaining open issue is why the atoms initially distributed in the whole lattice gather 
at a single site. Though it is clear that motion takes place via a tunneling process of atoms 
across the lattice, the way the atoms tunnel may involve not only single-order tunneling 
(tunneling of a single particle) but also higher orders of tunneling (tunneling of bound 
particles) j^ . 




Up to now, the analysis has been focused on the manipulation of atoms at the level of 
Floquet states. Clearly, to steer atoms in the lattice, one needs first to populate those states. 
Population of the Floquet state where atoms appear equally distributed across the lattice 
is realized by "slowly" turning on the periodic field, using a uniform distribution of atoms 
as the initial condition. That is, we ramp f(t) amplitude (Eqj5J), h, from to its maximum 
value. This process is carried out at a speed Of that satisfies the constraint af/u <C 1. 

In the following, we consider a lattice of 20 sites. In Fig]5]the atom number distribution is 
recorded after ramping up go during the two consecutive stages, as depicted in FigJU Notice 
the apparent localization of atoms on the sites that are perturbed by the AC pulse U(go) 
[see Eqj3]. In the first stage, the AC pulse is applied on the sites 2 and 3 and localization 
takes place on the site 3. In the second stage, the AC pulse acting on sites 13 and 14 causes 
localization on site 13. The localization on the perturbed sites depends on 9 similarly as in 
FigH] for the two modes system. 

A question that arises is whether this procedure can be applied to much larger lattices, 
and which limitations exist. Further tests indeed indicate that the movement of atoms is 
possible for a lattice with a larger number of sites, at the expense of reducing the ramping 
speed. The larger the lattice becomes, the longer the time is needed for the atoms to move 
across the lattice and gather together at few sites. 

From the quasienergy analysis, more sites in the lattice implies more Floquet states in 
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FIG. 5. (Color online) Atoms density \ipi\ 2 vs. location site i. (Upper) Q = 0; (Middle) Q = 0.025; 
(Lower) Q = 0.05; The dashed lines indicate the impurities location at the sites i = 8, 17. (see 
text for details on the localization sites) C = 1, oj = 4.7T, h = 4, af = 0.1, 9 = tt, a = 5 x 10~ 3 , 
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5. Computations are performed with the U(go) shape depicted in the upper inset of Figj3j 
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the quasienergy spectrum with energies becoming more densely packed and with smaller 
energy gaps. This increases the probability for Landau-Zener transitions, thus decreasing 
the control over the atoms. 

We have showed for a two-mode system that localization may take place against a bias. 



Lattice impurities cause a simi 
setting (cf bifurcation in Ref. 



ar effect in the energy spectrum to that of a bias for a dimer 



27| with inset in Figj2]). So, it is of particular interest to see 



what happens with the movement of atoms from one site to any other site in the presence 
of impurities. 

Here, we consider two impurities placed between the initial and final sites of the motion 
process. The impurities are introduced by changing the energies at some specific sites, i.e, 
E n = Q * (5 n ^ + 5 n m ), where 5 n ^ is the Kronecker's delta function, k and m are the indices 
of the impurity sites in a lattice that runs from n = 0, 20 and Q is the amount of energy 
at those sites. For low Q, movement of atoms takes place as in the homogeneous lattice, 
whereas for large Q, the atoms tend to localize at the impurity sites, shown in FigHJ 

We performed similar computations with a random distribution of impurities (not shown). 
The results observed are similar to those exhibited above, where the atoms tend to gather 
at impurities with large E n . 



C. Conclusions 



In summary, we have shown the steering of a Bose-Einstein condensate with a generalized 
management of the atomic interaction in a nonlinear discrete AC driven lattice. Transport 
of atoms is realized via the control of the atom's localization process. The whole process 
relies on gathering atoms in different sites of the lattice upon successive ramping move in 
time of a DC part of the atomic interaction strength. 

Localization of atoms is achieved upon ramping the atomic interaction assisted by an 
AC field, together with periodic oscillations of the nonlinear interaction. This process of 
localization is first explained within a two-mode approximation, where the sign of the popu- 
lation imbalances is controlled by the phase of the periodic functions. This is supported by a 
symmetry analysis. We also show within a two-mode setting that control of the localization, 
created by the application of AC perturbations, may overcome a bias. 

This was later shown to be feasible in the lattice, when movement of atoms was proven 
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to be robust against a low intensity of disorder. 

The control procedure exposed here paves the way for further application in other non- 
linear systems, such as nonlinear optics where the nonlinearity strength is modulated in 
optical waveguides. Likewise, these results can be used for a generic equation like the dis- 
crete self-trapping equation with application to fundamental problems, such as dynamics 



of small molecules 
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dynamics of molecular crystals [29]], amongst others. On the other 
hand, similar management of the atomic motion may be potentially useful in the control of 
momenta states in a condensate and consequently for the transport of ultra-cold atoms in 
optical lattices. 
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